Dynamic response of a mesoscopic capacitor in the presence of strong electron interactions 
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We consider a one dimensional mesoscopic capacitor in the presence of strong electron interactions and 
compute its admittance in order to probe the universal nature of the relaxation resistance. We use a combination 
of perturbation theory, renormalization group arguments, and quantum Monte Carlo calculation to treat the 
whole parameter range of dot-lead coupling. The relaxation resistance is universal even in the presence of 
strong Coulomb blockade when the interactions in the wire are sufficiently weak. We predict and observe a 
quantum phase transition to an incoherent regime for a Luttinger parameter K < 1/2. Results could be tested 
using a quantum dot coupled to an edge state in the fractional quantum Hall effect. 
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PACS numbers: 85.35.Gv, 73.21.La, 73.23.Hk, 73.43.Jn 

The dynamical response of mesoscopic conductors consti- 
tutes a mostly unexplored area of coherent quantum transport, 
which has recently led to groundbreaking experiments d. 
The mesoscopic capacitor \2] is one of its elementary build- 
ing blocks: a quantum dot influenced by an AC gate voltage, 
which is put in contact with an electron reservoir. It has been 
studied so far at the single electron level, with possible mean 
field generalizations fl. Both the capacitance and the re- 
laxation resistance R q , obtained from the low frequency ex- 
pansion of the admittance G(w) » -icoC^ + u) 2 C 2 R q , are fun- 
damentally affected by the quantum coherence of the device. 
At zero temperature, a single spin polarized channel yields a 
relaxation resistance R q = h/(2e 2 ), which is independent of 
the dot-reservoir connection. Ref. (H has confirmed this re- 
sult for a quantum dot with weak charging energy. 

However, quantum dots with reduced size exhibit strong 
Coulomb blockade, and there is also a clear need to analyze 
whether electron-electron interactions in the lead are relevant. 
Here, taking rigorously these aspects into account, we prove 
that there is quantum phase transition from a coherent to an 
incoherent regime, where a relaxation resistance cannot be 
defined. For weak interactions, the universal behavior is re- 
covered even in the presence of strong Coulomb blockade. 

We consider a quantum dot (Fig. [TJ connected to a reser- 
voir modelled by a Luttinger liquid lead, which allows to ac- 
count exactly for Coulomb blockade effects. We discuss sepa- 
rately the absence (Luttinger parameter K — 1) or the presence 
(K < 1) of interaction in the adjacent lead. This setup and the 
underlying physics is similar to that studied in Ref. j4j], where 
attention was solely focused on the static occupation of a res- 
onant level. Here we show that below K = 1/2 the Kosterlitz 
Thouless type phase transition driven by the dot-lead tunnel- 
ing strength triggers a transition of dynamical transport from 
a coherent to an incoherent regime, hence provoking a devia- 
tion from the universal R q = h/(2e 2 ). We use a combination 
of analytical (perturbation theory, renormalization group) and 
numerical (quantum Monte Carlo) approaches to monitor the 
capacitance and the relaxation resistance over the whole range 
of dot-lead connection. The present results can be applied to 



carbon nanotube quantum wires as well as dots defined in the 
fractional quantum Hall effect (FQHE). 

The starting point is the Hamiltonian for a non-chiral, semi- 
infinite Luttinger liquid where the dot region corresponds 
to the interval [0, L]: 
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The first part is the kinetic part, followed by the backscat- 
tering term at x — (strength V), and finally the contri- 
bution from the charging energy with Eq = e 2 /(2C) (C is 
the geometrical capacitance). The canonically conjugated 
fields </> and 9 satisfy the commutation relation [<p{x), 6{x')] = 
(/7r/2)sgn(x - x'). V is the backscattering strength on the 
point contact. Eq = e 2 /(2C) denotes the charging energy 
The time dependent gate voltage oscillates around V g . Us- 
ing the Matsubara imaginary time path integral formulation, 
the quadratic degrees of freedom away from x — can be 
integrated out. The kinetic part of the effective action then 
reads S kin = {nK[3Y l ^ |w„|/(l - e - 2rf ^ A )|0( W „)| 2 , where 
4>(u>„) is the Fourier transform of </>(t) (now specified at x — 0), 
and A = m'f/L is the level spacing. The same action can 
be derived alternatively starting from a single chiral Luttinger 
liquid "loop", hence the relevance for the FQHE regime Jg]. 
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FIG. 1: Left: schematic view of the mesoscopic capacitor: a ID 
quantum dot, capacitively coupled to a gate with time-dependent 
voltage V g + 6V g (t). Right: schematic phase diagram in the degener- 
ate case. V c denotes the critical backscattering strength for K = 1/3. 
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Within linear response in the oscillating gate voltage, the 
(imaginary frequency) admittance can be expressed as: 



G(ia) n ) = 



e 2 2\cj n 
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The dynamical conductance is obtained by analytic continua- 
tion G(a>) = G(iu>„ — > to + id), while the capacitance reads: 
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dr0(r) . (3) 



We start with a discussion of the weak barrier case, using 
perturbation theory in V/D (bandwidth D). The capacitance 
and relaxation resistance are derived as an expansion in orders 
of V, C, = cf + cf + Cf + ... and fi, = <> +*f + .... 
Introducing 
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one obtains to zeroth, first and second order: 



G '(io)„) = — 

h n pa„ 

„2 



G m {iu n ) = -t-^Vy/FJOj^ cos(2^iV) 



where we defined: 



r 

Jo 



F ± (v) = exp 
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dv [cos(47rA?)V + (w,„ v) + V-((o n , v)] (9) 



V ± (w„, v) = 1 - F ± (v) (1 ± cos(w„v)) (10) 
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The 2„ in F± is limited by D. From Eqs. ©-(IT]), the capaci- 
tance at low temperature becomes: 
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It is clear from these expressions that the total capacitance 
is a periodic function of N, with period 1 . Below we focus on 
the interval < N < 1 . The results for the relaxation resis- 
tance, at low temperature, are simple since the computation of 
the first and second order contribution shows that they vanish: 
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FIG. 2: (color online) Capacitance as a function of the backscat- 
tering strength V, obtained at temperature T = 0.04 with Monte 
Carlo computations (dashed lines). The solid lines show the pre- 
dictions of the perturbative calculations up to second order. We take 
Ec/tt 2 as the unit of energy and use the parameters D = 2nJ/p = Sn 



and A/(2K 2 ) 



Temperature dependence is shown in the bot- 



tom right panel, where the vertical axis is measured on a logarithmic 
scale. 



The charging energy thus does not modify the value of relax- 
ation resistance, while electron interactions in the lead intro- 
duce a trivial factor 1 /K . At zero temperature, the sums and 
integration of Eqs. (|5])-(|7| can be done analytically in certain 
cases. For example, when Ec = and K = 1, one has F + (Q) = 
(A/(2?rD)) 2 {D » A), and one can show that the result for the 
capacitance is = (e 2 /A)(l - 2rcos(2nN) + 2r 2 cos(4ttN)), 
with r - nV/D. This coincides with the development of the 
non-interacting formula found in Ref. 01 in powers of the re- 
flection coefficient r. = (e 2 /A)(l - r 2 )/(l - 2rcos(2nN) + 
r 2 ). In the more general case of non-zero Ec, and K + 1, one 
has F + (0) ~ (E t K/(ttD)) 2k , and the integration of Eq. Q has 
to be computed numerically. 

The perturbation theory thus proves the universality of the 
charge relaxation resistance in the weak barrier limit even in 
the presence of interactions. To study the non-perturbative 
regime, the path-integral Monte Carlo method is applied to the 
action for the discretized path </>(t = j/3/ J) (J = 0, 1, • • • , J - 
1). We estimate thermal average by generating discretized 
paths using local update in the Fourier space and the clus- 
ter update |7. ;&]. The top (bottom) row of Fig. [2] shows the 
calculated capacitance as a function of V at T = 0.04 and 
K = 1 (K = 1/3). The left and right columns correspond to 
the non-degenerate case (N = 0) and the charge degenerate 
case (N = 1 /2), respectively. With increasing V the Coulomb 
staircase becomes sharper, which results in the decrease (in- 
crease) in Cf, = d(Q)/dV g in the case of N = (N = 1 /2). 
The second-order perturbation theory, shown as solid lines, 
displays an excellent agreement for small V. Especially, it is 
remarkable that only for the case of K — 1/3 and N — 1/2 (the 
right bottom panel of Fig. [2J, exhibits an abrupt increase at 
a finite V, signaling a possible transition. One can see that C M 
grows as oc 1/T in the large barrier region. 
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To reveal the origin of the transition behavior, we next 
examine the strong barrier limit using an instanton method 
which was developed for the Kondo model J^t]. Near the de- 
generacy point N - 1/2, the configuration of the bosonic field 
(f> can be represented in the dilute instanton gas approximation 
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where Sj = s(—iy , and s ~ 1 denotes the separation be- 
tween the well minima (0 is the step function). Inserting Eq. 
(TToT l in the full effective action, the partition function becomes: 

Z = y t 2 " dr 2 „ dr^ • • • dn 
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where t is the tunneling amplitude between the well minima, 
and t c is the short-time cutoff, u - {2N - \)Ej denotes the 
deviation from the degeneracy point. Note the similarity be- 
tween this partition function and that which was proposed in 
the context of dissipative Josephson junctions [| 1 Ofl - One can 
therefore identify the scaling equations [12]: 
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which are familiar in the context of a Kosterlitz Thouless 
transition in the two-dimensional XY model. No further ar- 
guments are needed when one deviates from the degeneracy 
point: since t is small, starting from u + 0, Eq. < fT~9b predicts 
that u will further increase, leading the system further away 
from the degeneracy point. This means that <p will be trapped 
in an effective harmonic potential, and one thus recover the re- 
sult of Eq. (0, which is therefore universal. For the charge de- 
generate case N - 1/2, the transition corresponds to a Kondo 
type transition associated with the charge pseudo spin on the 
dot. Eqs. ([T8l determine the tendency of the dot-lead trans- 
mission as temperature is lowered; (f, s 2 ) flows along one of 
the hyperbolic curves B 2 -4t 2 = const., where B = l-s 2 /(2K). 
For K > 1/2, the tunneling strength always grows upon re- 
ducing the temperature, and the system reaches the Kondo 
fixed point where the dot is strongly coupled to the reservoir. 
An electron freely tunnels in and out of the dot irrespective 
of the initial tunneling strength. In particular at K = 1 this 
implies that the charge relaxation resistance is universal, i.e., 
R q = h/(2e 2 ), as a consequence of the unitary limit of the un- 
derlying Kondo model. On the other hand for K < 1/2, there 
is the possibility that at a critical, sufficiently weak transmis- 
sion t ("large" V), the RG flow always drives the system into a 
weak coupling configuration with specified charge. Then the 
charge fluctuation remains finite, i.e., Q> 2 ) - {(j)) 2 - (ns/2) 2 , 
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FIG. 3: (color online) Extrapolated value R(ko„ — > 0) as a function 
of V at temperature T = 0.04 in the degenerate case. For D and 
A/(2K 2 ), we use the same parameters as in Fig. [2] Inset: R(ia>„ — > 0) 
for K = 1/3 at T = 0.01(o), 0.04(A) and 0.08(o) in the vicinity of the 
location where the crossing occurs. 



so that the capacitance diverges as oc l/T at low temperatures 
[see Eq. {3)]. This explains the transition observed for the 
capacitance in the strongly interacting case. 

We now describe the effect of the KT transition on the dy- 
namical properties. If u> <K 1/trc holds (with RC time trc), 
the charge relaxation resistance can be defined in the low- 
frequency expansion G(w) = -iuC^ + u 2 C 2 Rq + 0(ca 3 ). How- 
ever, the validity of this expansion is not obvious, since the KT 
transition may influence trq itself. Instead, we investigate the 
low-frequency resistance using 



R(ia)„) = 
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where G(ia> n ) and are defined in Eqs. (f2]i and (f3), respec- 
tively. The extrapolation R(ia>„ — > 0) gives the real part of the 
impedance in the low-frequency limit, hence trc = R(0)Cft. 
In Fig. [3] we plot R(0) for K — 1,1/3 and 1/5 as a function 
of V at temperature T = 0.04. For K = 1, R(0) equals h/(2e 2 ) 
irrespective of V, in agreement with the universal charge re- 
laxation resistance |2|,|3j]. For K — 1/3 and 1/5, the univer- 
sality is observed in the weak barrier region, whereas R(0) is 
abruptly enhanced with increasing V, reflecting the RG flow 
to the weak coupling regime due to the KT transition. The 
temperature dependence of R(0) for K = 1/3 is shown in the 
inset of Fig [3] which indicates that R(0) diverges as T — > in 
the strong barrier region. 

The KT transition plays a crucial role in the relevance of 
the universal charge relaxation resistance. If 2t + B > 0, the 
system scales to the weak barrier limit, where trc is indepen- 
dent of temperature. If 2t + B < 0, on the other hand, the 
scaling equations ( TT8l predict s 2 — > const, and t oc T~ B , so 
that trc roughly scales as oc j-\(j 2B + const.), which grows 
faster than the (thermal) coherence time r co h oc 1 /T as temper- 
ature is lowered. These observations suggest that if 2t + B > 
coherent transport can be realized by lowering temperature to 
guarantee trc < T" co h, while if 2t + B < electronic transport 
in the dot decoheres before charge relaxation is achieved. In 
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FIG. 4: (color online) Product G{icQ„ = 0)R(ico n = 0) for K = 1/3 as 
a function of the backscattering strength V, for different temperatures 
T. For D and A/ (2K 2 ), we use the same parameters as in Fig.[2j Inset: 
^(OJC^r as a function of V, for different T; the crossing of the curves 
gives a good estimate of the transition point. 



the latter case, the quantum dot effectively acts as a reservoir 
and consequently the dynamical property of the system is gov- 
erned by transport through the point contact between the two 
"reservoirs". Therefore the V-dependent low-frequency resis- 
tance observed in the inset of Fig.[3]reflects the revival of the 
Landauer-type transport. To see this behavior more clearly, 
we plot in Fig. [4] the product G(ico n — > Q)R(ia>„ — > 0) for 
K — 1/3 as a function of V. In the strong barrier region, 
G(0) is finite and equal to [/?(0)]~', which is a familiar prop- 
erty of transport through a point contact. Upon decreasing V, 
however, G(Q)R(Q) is suppressed since G(0) decays to zero be- 
cause of charging up, although R(0) — > R q is finite. Moreover, 
we see that the coherent region G(0)R(0) = extends to larger 
V upon lowering temperature. Finally, we determine the phase 
boundary of the coherent-incoherent transition by tracing the 
temperature dependence of the ratio T^c/Tcoh = RifyC^T. The 
above discussion suggests that there exists a critical backscat- 
tering strength V c , below which T/jc/fcoh decays to zero, while 
it diverges otherwise (see the right panel in Fig. [TJ. From the 
inset of Fig. 21 the critical value is estimated as V c - 7. 

In conclusion, the study of the mesoscopic capacitor in the 
presence of strong electron electron interaction shows that the 
relaxation resistance for a dot connected to Luttinger liquid 
is universal R q = h/(2e 2 K) as long as interactions are suffi- 
ciently weak. Below K < 1/2, this resistance is governed by 
the strength of the dot-lead coupling: at the charge degener- 
acy point, there is a critical coupling strength, governed by a 
KT type phase transition, below which the dot acts as an inco- 
herent reservoir and the low-frequency resistance exceeds the 
universal value. In this incoherent regime, the charge relax- 



ation resistance cannot be defined anymore due to the diver- 
gence of the RC time. 

Results could be probed experimentally using quantum dots 
connected to an edge state in the FQHE regime. Another ex- 
perimental probe could use one dimensional quantum wires 
(non chiral Luttinger liquids) with the limitation that the oper- 
ating frequency would have to be larger than the inverse time 
of flight within the wire, in order to avoid renormalization 
effects due to eventual Fermi liquid leads connected to this 
wire lllll . 
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